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ABSTRACT 

We present spectroscopy of 28 SNR candidates as well as one H II region in M81, and two SNR 
candidates in M82. Twenty six out of the M81 candidates turn out to be genuine SNRs, and two in M82 
may be shocked condensations in the galactic outflow or SNRs. The distribution of [N Ilj/Hct ratios 
of M81 SNRs is bimodal. M81 SNRs are divided into two groups in the spectral line ratio diagrams: 
an [O III]-strong group and an [O III]-weak group. The latter have larger sizes, and may have 
faster shock velocity. [N II]/Ha ratios of the SNRs show a strong correlation with [S II]/Ha ratios. 
They show a clear radial gradient in [N II]/Ha and [S II]/Ha ratios: dLog ([N II]/Ha)/dLog R 
= —0.018 ± 0.008 dex kpc“^ and dLog ([S II]/Ha)/dLog R = —0.016 ± 0.008 dex kpc“^ where R 
is a deprojected galactocentric distance. We estimate the nitrogen and oxygen abundance of the 
SNRs from the comparison with shock-ionization models. We obtain a value for the nitrogen radial 
gradient, dLog(N/H)/dLogR = —0.023 ± 0.009 dex kpc“^, and little evidence for the gradient in 
oxygen. This nitrogen abundance shows a few times flatter gradient than those of the planetary 
nebulae and H II regions. We find that five SNRs are matched with X-ray sources. Their X-ray 
hardness colors are consistent with thermal SNRs. 

Subject headings: galaxies: spiral — galaxies: individual (M81, M82) — galaxies: ISM— galaxies: 
abundances 


1. INTRODUCTION 


Supernova remnants (SNRs) in nearby galaxies are 
an ideal target to study various statistical aspects of 
SNRs, because we can find them in the entire region of a 
galaxy (within the survey limits), and because they are 
in the same distance from us. SNR candidates in nearby 
galaxies are often found using narrow band imaging with 
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erences therein). Optical spectroscopy of the SNR can¬ 
didates is a powerful tool to confirm whether they are 
genuine SNRs or not, to derive their physical parame¬ 
ters (electron densities, temperatures, shock velocities) 
and chemical abundances, an d to investigat e the radial 
gradient of the abundances (IDoo ita 1 1197 ^ [Raymond I 
119791 : IShull h McKeel[T979l: IDooita et al.lll984h . OoticM 
spectroscopic studies of SNRs in several nearby galaxies 
have been done in the past. The Local Group galaxies 
have been primary targets of optical spectroscopic stud¬ 
ies: LMC, SMC (iRn ssell fc Dooital 119901: iPa vne et al.l 
[200l), M31 (IBlair et al.lll981L 1 19821:[Galarza et al.lll999n 
and M33 (iBlair fc Kirschner 71985 : ISmith et al 119931 : 
iGordon et al.lll998f) . Recently SNRs in several low lumi- 
nosity disk galaxies ha ve been studied spectroscopically 
(iLeoiridaki et al.ll2013^ . 

IGordon et al.l ()1998[1 performed an extensive optical 
study for a large sample of SNRs in M33 includini 
those in the previous studies (|Blair fc Kirschne^ 119851 
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ISmith et al.l 1199311 . They presented several interesting 
results on the correlation between line ratios: no relation 
between [N II]/Ha and [S II]A6717/[S II]A6731, a weak 
relation between the SNR diameter and [N II]/Ha as 
well as [S II]/Ha, and a strong correlation between 
[N II]/Ha and [S II]/Ha, and between [O I]/Ha and 
[S Il]/Ha. They found also that the line ratios of 
[N II]/Ha of M33 SNRs decrease as the galactocentric 
distance (GCD, R) increases at 0 < i? < 6 kpc, and they 
show a significant dispersion at a given galactocentric 
distance. They suggested that the significant dispersion 
is mainly due to a larger abundance dispersion in the 
inner region of the galaxy. However, their spectra cov¬ 
ered only the red region, 6200 to 7500 A so that they 
could not have any information of the blue emission lines, 
exi Optical sp ectroscopic s tudies of th e SNRs in M31 
were given by IBlair et HI (|198lL 1198211 : iGalarza et al.l 
(1199911 . From the analysis of the co mbined sample in¬ 
cludi n g those in the previo us studies (IBlair et al.l I1981L 
1198211 . IGalarza et al.l (1199911 found a clear radial gradient 
of [N II]/Ha in the range of 3 < i? < 20 kpc (—0.04±0.01 
dex kpc“^ in their Fig. 6d), and little evidence for gradi¬ 
ents in [S II]/Ha, [O III]/H/3, [O III]/[0 II], and electron 
temperature (Tg). They found also a significant corre¬ 
lation between [O III]/H/3 and [O II]/H/?, noting that it 
may be mainly due to oxygen abundance variation. 

M81 (NGC 3031) is another nearby bright spiral 
(SA(s)ab) galaxy, located beyond the Local Group. It 
is a main member of the M81 Group, interacting with 
the well-known starburst galaxy M82 an d a dwarf galaxy 
NGG 3077 with morphology type of 10. iMatonick et al.l 
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(|1997fl provided a list of 41 SNR candidates selected us¬ 
ing the criterion of [S II]/Ha > 0.45 from the [S II] and 
Ha images of M81. However, they obtained optical spec¬ 
tra of only four in their sample (MF No. 2, 17, 18, and 
25) covering 4800-7200 A, and conhrmed that they are 
indeed SNRs. 

In this study we present a spectroscopic study of a large 
number of SNR candidates in M81, using a wide spec¬ 
tral range to cover from [O II]3727 to [S II]A6731 lines. 
This paper is composed as follows. Section 2 describes 
how to select and observe the targets and how to reduce 
the spectroscopic data. In §3 we classify SNRs, inves¬ 
tigate the property of emission line ratios of the SNRs, 
and compare them with shock-ionization models. Then 
we investigate the relations between the line ratios and 
sizes of the SNRs, and any radial variation of the line ra¬ 
tios and abundances. In §4 we discuss the radial gradient 
of chemical abundances derived from the SNRs in com¬ 
parison with those based on H II regions and planetary 
nebulae (PNe) in M81. We compare the radial gradients 
of M81 SNRs with those of the SNRs in other nearby 
galaxies. Also we match optical SNRs with X-ray source 
catalogs of M81. Final section summarizes the main re¬ 
sults with conclusions. 

We adopted a d istance to M81, 3.63 ± 0.14 Mpc 
(|Durrell et al.l 1201011 . inclination angle, z = 58 deg, 
and position angle, PA = 157 deg to calculate the 
depro iected galactocentric dista, n ces (R) of the tar - 
gets i|de Vauconlenrs et al.l 119911 : iPuerari et al.l 120141 1. 
The standard radius {R 25 ) of M81 is 18.34 kpc 
and i ts heliocentric radial velo city is — 34 ± 4 km 
s“^ (Ide Vauconlenrs et al.l 1199111 . One arcsec corre¬ 
sponds to 17.6 pc at the dist ance of M81. Th e distances 
to M81 and M82 are similar (|Lim et al.ll2013ll . 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Target Selection 

Primary targets (20) were sele cted from the lis t of th e 
SNR candidates in M81 given bv iMatonick et all (|1997[1 . 
We noted their IDs starting with MF. Then we added 
eight secondary targets we selected as SNR candidates 
from the Ha and [S II]AA6716,31 images obtained with 
INT 2.5m in the Isaac Newton Group Archive. Their IDs 
run from LI to L8. We also selected one H II region in 
the southeast region of M81 for reference. 

In addition, we checked the Ha and [S II] images 
of M82 in the Hubble Space Telescope (HST) archive 
(PI: Mountain and PID: 10776 for Ha images, and 
PI: O’Connel l and PIP: 11360 for [S II] images). 
Ide Criis et al.l ()2000tl presented a list of 10 SNR candi¬ 
dates in the northeast starburst region (M82B) of M82 
that are compact Ha -bright sources, selected from the 
Ha images obtained with HST/WFPC2. There are no 
spectroscopic studies of these candidates so that their 
nature remain to be known. Similarly we selected two 
compact Ha -bright objects located in the outflow re¬ 
gion of M82 as interesting targets. Figures [T] and 
display the location of the targets in M81 and M82, re¬ 
spectively. Note that the M81 targets are mostly located 
along the spiral arms. 

2.2. Observation 
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Fig. 1.— The spatial distribution of the SNR candidates (circles) 
and one H II region (diamond) in M81 in this study, marked in the 
grayscale map of the digital sky survey. North is up and east to left. 
Filled circles represent [O IIl]-weak SNRs (with [O IIl]/ H/9 < 1.0), 
and filled squares SNRs matched with X-ray sources IlSell et al.l 
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Fig. 2.— Location of two SNR candidates in M82 (M82L1 and 
M82L2) in this study, overlayed in the gray scale map of continuum- 
subtracted ACS Ho; image in the HST archive (PI : Mountain, PID 
: 10776). North is up and east to left. Note that they are located 
in the outflow direction. 


To obtain the spectra of the targets we used the Hec- 
tospec fiber-fed spectrograph ( Hectospec) on the 6.5m 
Multi Mirror Telescope (MMT) (|Fabricant et al.l [20051) . 
We selected a 270 mm“^ grating at a dispersion of 1.2 
A pixeU^, covering 3650 A to 9200 A at a spectral 
resolution of about 5 A. However, the sensitivity at the 
blue and red limits are low so that we used only the 
spectral range of 3700 A to 7000 A for analysis. The 
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spectrograph feeds 300 fibers over one degree diameter 
field of view and the diameter of each fiber is ^ 1"5, 
corresponding to 26.4 pc at the distance to M81. 

The observations of the emission line objects in M81 
and M82 were carried out as a part of our larger M81 
program in queue mode during the months of February 
to May 2014. We took 3 to 5 exposures of 1440 s for 
one configuration. We observed four configurations that 
are overlapped partially with each other. Total exposure 
times for one configuration are 4320 s to 7200 s. Seeing 
values during the observations ranged from I'.'O to 1'.'2. 


2.3. Data Reduction 

We reduced t he data using HSRED version 2 writ¬ 
ten by R.CooQ (|Kochanek et al.1 [20121) . including debi¬ 
asing, flat-fielding, aperture extraction of spectra, wave¬ 
length calibration, and flux calibration. Sky subtraction 
for spectra extraction was done using the average of the 
spectra for the nearest blank-sky fibers. 

We measured the emission-line fluxes using Gaussian 
fitting after continuum subtraction in IDL. The errors for 
the fluxes were estimated including the systematic errors 
of the flux calibration, background estimation, and sky 
subtraction. For extinction correction we derived c(H/l), 
the logarithmic extinction at H/1 , adopting the intrinsic 
value of Hct /H/3 = 3.0 for the SNRs based on the shock 
models (jRaymond II1979I : IShull fc McKedllQTfK l and 2.86 
for HII regions based on Case B recombination a nd a 
gas with Te r::; 10^ K (jOsterbrock fc Ferland II2006D. We 
adopted the extinction curve bv iCardelli et al.l ( 198911 . 

We measured the radial velocities of the targets us¬ 
ing HSRED version 2, after checking visually whether 
each measurement is reasonable or not. We adopted the 
range of 3900A to 6400 A for velocity measurements. 
The mean value of the measurement errors of the radial 
velocities is 3 km s“^. Table [T] lists a list of spectroscopic 
parameters of the targets: ID, R [kpc], radial velocity 
with its error, c(H/3), H/3 strength, and the extinction 
corrected line fluxes and errors, which are normalized by 
H/3 = 100. Hereafter we use only the extinction-corrected 
values for the line fluxes and line ratios for the following 
analysis. 

One of the secondary targets (L7) showed too poor a 
spectrum to use. Another of the secondary targets (L8) 
turned out to be a background galaxy at z = 0.33. Thus 
we used the spectra of 29 objects for analysis. Two of 
our sample are matched with the SNRs for whi ch spec¬ 
troscopic data are given in iMatonick et al.l (jl997f l: MF17 
and MF25 . The line ratios of t hese two SNRs in this 
study and iMatonick et al.l (1199711 agree well within the 
differences of about 20%. 


3. RESULTS 
3.1. SNR Classification 

Figure [3] displays stamp maps of the SNR candidates 
and one H II region in M81, based on the continuum- 
subtracted Ha and [S II] images we prepared using the 
INT images. We estimated visually the diameters of the 
SNR candidates in the continuum-subtracted [S II] im¬ 
ages, as marked by circles in EigurelU and listed them 
in Table O Note that the contrast between the Ha and 



Fig. 3.— Stamp maps (20^^ x 20^^) of the SNRs and one H II re¬ 
gion in M81 in this study, made from the continuum-subtracted 
Ha and [S lijimages. North is up and east to the left. Circles 
represent the SNR size estimated in this study. Note that the 
fiber with diameter 1^5 was located in the center of each map. 
Five SNRs are matched with the X-ray source catalog (ISell et al.l 
[20TJ): MF10=Sell259, MF17=Selll93, MF19=Selll95=Liul578, 
MF25=Selll72, and L5=Sell50=Liul653. 
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Fig. 3.— Continued. 
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Fig. 3.— Continued. 



Fig. 4.— Same as Figure [3] but for M82 SNRs. Each map cov¬ 
ers 6^^ X 6^^ made from the continuum-subtracted ACS Ho and 
[S II] images. Circles represent the size of the SNRs. Note that 
they show shell structures typical for SNRs. 


[S II] images is much higher for the H II region than for 
the SNR candidates. In Figure|^we showed stamp maps 
of two M82 SNR candidates, based on the continuum- 
subtracted Ha and [S II] images we prepared using the 
ACS Ha and [S II] images in the HST archive. Note that 
the high resolution images of these two objects show shell 
structures typical for SNRs or shocked regions. 

Figures [S] and [6] show spectra of 26 M81 SNR candi¬ 
dates, two M82 SNR candidates and one M81 H II re¬ 
gion, respectively. Note several prominent lines typi¬ 
cal for SNRs: [O II]AA3727,9, H/3 , [O III]AA4959,5007, 


[O I]A6300, [N II]AA6548,83, Ha , and [S II]AA6717,31. 
[O III]A4363 lines are too weak to measure for all objects. 
Hereafter we use for common doublet lines, [O II] for 
[O II]AA3726,9, [O III] for [O III]AA4959,5007, [N II] for 
[N II]AA6548,83, and [S II] for [S II]AA6716,31. 

The [S II]/Ha ratio has been considered to be an ex¬ 
cellent cnterion to distinguish SNRs from H II regions 
(|Long| |1985D. We derived the values of the JS II]/Ha ra¬ 
tio of M81 SNR candidates using their fluxes measured 
in the continuum-subtracted Ha and [S II] images. In 
Figure (3 we compared thus photometrically derived 
values with those we measured from the spectra. It 
shows a reasonable correlation between the two values, 
but with a large scatter. Linear-least-squares fitting 
yields [S II]/Ha (spec) = (0.894±0.059) (S II]/Ha (phot) 
(0.121 ± 0.061) with rms = 0.084, showing that its slope 
is slightly smaller than one. This confirms that the pho¬ 
tometric method based on [S II] and Ha images to search 
for SNR candidates in nearby galaxies is an efficient tool. 

We plotted the distributions of the [S II]/Ha and 
]N II]/Ha ratios of M81 SNR candidates derived from 
spectra in Figure El We a lso plotted the data for M31 
SNRs (jGalarza et ffi.l [19^ . M33 SNRs (jGordon et al.l 
FT^. a,Tid SNH s in six low luminosity disk galaxies 
(|Leonidaki et al.l 1201311 for comparison. These figures 
show two notable features. First, the [S II]/Ha ra¬ 
tios of all 26 SNR candidates in M81 are higher than 
0.5, showing that they are indeed SNRs. Second, the 
[N II]/Ha distribution of M81 SNRs is bimodal with 
peaks at 0.75 and 1.15. M31 SNRs also show a simi¬ 
lar bimodal distribution. However, the number ratio of 
the higher and lower components of M81 SNRs is larger 
than that for M31 SNRs. It is noted that the presence of 
the higher component is mainly due to the SNRs located 
in the inner region at i? < 5 kpc, as shown later. 

In Figure llKa), we displayed the radial velocity with 
respect to the galaxy center of M81 SNRs as a func¬ 
tion oHgositffin angle, in comparison with those for HI 
gas (lRotslll975ll and globular clusters (|Nantais fc Huchral 
l201(iL We derived the rotation velocity of the SNRs, 
using Vc = (vr — Vo)/ sinicos(^ — 4>o) where </> and 4 >o 
are position angles of the SNRs and the minor axis of 
M81, respectively, and i is an inclination angle. We 
ado pted the systemic velocity of M81, vq = —34 km 
s~^ (lde Vaucouleurs et al.iri991h . We plotted the results 
for the SNRs with {(f) — (f)o) > 20 deg as well as those of 
HI gas (jRotsl[T97^ in Figure (HFb). It is found that the 
rotational velocity data of M81 SNRs are in excellent 
agreement with those of HI gas, and are consistent with 
those of the globular clusters. This shows that all M81 
SNRs in this study are located in the disk, following the 
disk rotation, so that they belong to the disk population 
of M81. 


3.2. Emission Line Diagnostics 


General characteristics of emission lines typical for 
shocked regions such as SNRs are des cribed in rela¬ 
tion with s h ocked -i onization models jDogitaJJl^TTl 


Dooita et al.1 119841 IRavmondl 
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Shull V McKef 

1979|:IBlair et al.lll982llSmith et a 
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Matonick et al. 

19971 iGordon et al.l 119981: 

Galarza et al.l 119991 

Douita & Sutherland! I1995L 119961; lAllen et alJ I2008D. 


We give only a brief summary of them as follows. 
[O III]A5007 lines are emitted in the region close to 
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Fig. 5.— Spectra of M81 SNRs in this study. Note several prominent lines typical for SNRs: [O II]AA3727,9, H/3, [O III]AA4959,5007, 
[O I]A6300, [N II]A6548,83, Ha, and [S II]AA6717,31. 
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Fig. 6.— Spectra of one H II region in M81 and two SNR candi¬ 
dates in M82. 


the shock front so that their brightness is affected by 
the postshock electron temperature and the oxygen 
abundance. The line ratio of [O III]A5007/[0 II]A3727 
depends mainly on the postshock condition rather 
than on metallicity (iDopita I 119771) . The line ratio 
of [O III]AA4959,5007/[O III]A4363 is an indicator of 


o 

u 

a, 

t/3 


t-- 

vo 


1.8 

1.6 

1.4 

1.2 

1.0 

0.8 

0.6 

0.4 



0.5 1.0 1.5 

[S II]/Ha (Phot) 


Fig. 7.— Comparison of [S Il]/Ha ratios of M81 SNRs derived 
from photometric images and spectroscopic data in this study. The 
dashed-line and dotted-line represent one-to-one relation and the 
best linear-fitting relation, respectively. 


electron temperature. [O III]A4363 lines in SNRs are 
weak in general, although they are stronger than in 
H II regions. [N II] lines are emitted in the large 
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Fig. 8.— The [S Il]/Ha (a) and [N Il]/Ha (b) ratio dis- 
tributi ons of M81 SNRs (this study) i n comparison with M33 
SNRs IIGordon et al.iri998li, M31 SNR s IjGal arza et aTIllOOOli and 
SNRs in six disk galaxies IILeonidaki et al.ll2013l) . Note that the 
[N Il]/Ho ratio distribution of M81 SNRs is bimodal. 
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Fig. 9. — Radial velocity with respect to the galaxy center of 
M81 SNRs (circles) vs. Position Angle (a), and corresponding 
rotational velocity versus galactocentric distance of the SNRs 
with position angles larger th an 20 deg (b ). Triangles and crosses 
repre sent the data for HI gas fRotsHlOTSh and globular clusters in 
M81 IINantais fc Huchrall20inl ~l. Note that M81 SNRs follow well 
the rotation of the M81 disk. 


recombination region behind the shock front, and are 
known to be an outstanding tracer of nitrogen abun¬ 
dance, little affected by shock temperature or electron 
density. Thus the line ratio of [N II]/Ha is a good indi- 
cator of relative abundances of n i trogen and hydrogen 
(iDooita et al.l iSrnith et al.1 1 199,11 iGordon et al.l 

ll998HGalarza et al.lll999ll . The ratio of the doublet lines, 
[S II]A6717/[S I I]A6731, is a well-know n indicator of 
electron density (|Blair fc Kirschn^ll985D . The ratio of 
[S II]A6731/Ha is mainly sensitive t o abundance, being 
little affected by shock conditions ( Dopi ta et al.l 11984 
iBlair fc Kirschneilll98!Tl iSmith et al.l 1199, lib [O I]A6300 
lines become stronger as the shock velocity increases 
l|Russell fc Dopital 119901) . The ratios of these lines with 
respect to H,5 vary significantly depending on shock 
velocity in the low shock velocity range (if the shock 


Fig. 10.— The [O IIl]/H/3 vs. [O Il]/[0 III] diagram (a) and the 
[N II] /Ha vs. [O II] /[O hi] diagram (la) of M81 SNRs (open cir¬ 
cles), H II region (asterisk), and M82 SNR candidat es (crosses) in 
this s tudy in comparison with M31 SNRs /diamonds. IGalarza et all 
1199911 . Open and filled circles denote (O Ili]-strong and [O IIl]-weak 
SNRs, respectively. The dashed and solid lines represent the upper 
limit of the theoreti cal models of H II reg ions for one and two solar 
abundance given by IDopita et al.1 II2013I1 . 

velocity is smaller than ^ 100 km s“^). However, some 
of them ([O II]AA3727,9, [O III]A5007, [N II]A6584, 
and [S II]A6731) change little when the shock velocity 
is larger tha n ~_100 km s“^(see Figures 5 and 6 in 
iDopita et al.l (11984 )1 so that they can be used as an 
abundance indicator. 

Figure llOK al displays the relation between 
[O III]A5007/H/3 and [O II]/[O III]A5007 of the 

SNRs in M81 (this study) in comparison with those in 
M31 (jGalarza et al.l 119991) . We plotted also lines for 
the theoretical models for H II regions w i th on e and 
two solar abundance given bv iDopita et al.l (I2013D . The 
following features are noted in this figure. First, M81 
SNRs are located clearly in two groups: an [O III]-weak 
group ([O III]A5007/H/3 < 1) and an [O III]-strong group 
([O III]A5007/H/3 > 1). A majority of M81 SNRs (20, 
77 %) belong to the [O III]-strong group, and 23% (6) 
to the [O II I]-weak group (LI, L2, L4, MF 22, MF29 
and MF35). iDopita et al.l (11984 presented the strength 
of some emission lines as a function of shock velocity 
(vs) in the shock-ionization model for solar abundance 
in their Figures 5 and 6. According to these models, the 
value of the [O III]A5007/H/3 ratio changes rapidly from 
-0.5 at Vs ~ 80 km s“^ to -)-0.5 at Vs Ri 100 km s“^, and 
increases slowly to ~ 0.65 at Vs ~ 170 km s“^. Thus the 
[O III]A5007/H,5 ratio becomes degenerate for Vs < 80 
km s“^. On the other hand, the [O II]3737/H,5 increases 
more slowly from 0.25 at Vs = 50 km s“^, to 1.1 at 
Vs = 90 km s“^ so that it is a good indicator for shock 
velocity in the low shock velocity range. The measured 
line ratios of the [O III]-weak SNRs indicate that their 
shock velocity may be 65 to 80 km s“^. Second, the 
[O III]A5007/H/3 ratio shows a strong correlation with 
the [O II]A3727/[0 III]A5007 ratio, decreasing as the 
[O II]A3727/[0 III]A5007 ratio increases. M81 SNRs 
and M31 SNRs are located along the same sequence. 
However, the fraction of the [O III]-strong SNRs in M81 
is larger than in M31. This indicates that the metallicity 
of M81 SNRs as well as M31 SNRs may be higher than 
the low metallicity adopted for the H II region model. 
Third, two M82 objects belong to the [O III]-weak 
group. 

However, this is not the end of the story for the [O III]- 
weak SNRs. According to the modern shock models 
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Fig. 11.— The [N II] /Ha vs. [S II] /Ha diagram (a) and the 
[O l]/Ha vs. [S H]/Ha diagram (b) of M81 SNRs (circles), H H re¬ 
gion (asterisk), and M82 SNR candidates ( crosses) in this study , 
in comparison with M31 SNRs (diamond s, IGalarza et al.iri999f) . 
M33 SNRs (triangles, IGordon et al.1 [19981') . and SNRs in six disk 
galaxies (plusses, ILeonidaki et aI.1l2Q13l) . Open and filled circles 
denote [O Hl]-strong and [O Hl]-weak SNRs, respectively. Two 
dashed lines for [S II]/Ha =0.45 and [N II]/Ha =0.5 were used 
for selecting M31 SNRs in IGalarza et aO 1)199911 . and another for 
[N Ilj/Ha =0.63 is a lower limit for shock ionization models. In 
(a) the solid line represents a linear fit for M81 SNRs. The dotted- 
line represents a linear fit result for M81 SNRs excluding one point 
(at [S H]/Ha ~ 1.1 and [N H]/Ha 0.78). In (b) the solid line 
and dashed line represent linear fits for M8i SNRs and M33 SNRs, 
respectively. 


such as those in iHo et alJ (I2014D . the fast shocks would 
be [O III]-weak objects rather than [O III]-strong ob¬ 
jects, and they are expected to have strong [O I]/Ha, 
[N II]/Hab, and [S II]/Ha ratios. However, the [O III]- 
weak objects in Figures 11 and 13 show a large range of 
[O I]/Ha, [N II]/Hab, and [S II]/Ha ratios. This suggests 
that these objects may be in part photo-ionized by stars, 
and are in part ionized by fast shocks with 200-300 km 
s“^. 

Similarly we plot the [N II]/Ha versus 
[O II] /J O III]A5007 diagram of the SNRs in Fig¬ 
ure llOlf b'). It is found that the [N II]/Ha ratio shows 
no correlation with the [O II]/[O III]A5007 ratio. This 
indicates that the [N II]/Ha ratio depends little on shock 
conditions, in contrast to the [O III]A5007 ratio. Also 
M81 SNRs have, on average, higher [N II]/Ha ratios 
than M31 SNRs. This indicates that M81 SNRs may 
have higher nitrogen abundance compared with M31 
SNRs. Two M82 SNRs are located in the region of M31 
SNRs. On the other hand, one M81 H II region is far 
below the locations of the SNRs. 

Figure llllf al displays the relation between 
[N II]/Ha and [S II]/Ha ratios of the SNRs in M81 
(this study) a s well as those in M31, M33, and other 
disk galaxies (IGalarza et all Il999t IGordon et al.l 119981 : 
ILeonidaki et al.l 1201311 . We plotted also the boundary 
lines at [S II]/Ha =0.45 and [N II] /Ha =0.5, whi c h wer e 
used for selecting M31 SNRs by ICalarza et al.l (|1999H . 
and a line at [N II]/Ha = 0.63 (Log [N II]/Ha = —0.2), 
a lower limit for shock ionization models. Several 
features are noted in this figure. First, all M81 SNRs 
are located above the lower limit for shock ionization 
models, while some of the SNRs in other galaxies 
are below this limit. Second, M81 SNRs show a 
remarkably strong correlation between these ratios, 
which is fitted linearly well with a slope close to one: 
[N II]/Ha = (0.934 ± 0.080) [S II]/Ha -b(0.112 ± 0.071) 
with rms = 0.109 ([N II]/Ha = (0.979 ± 0.066) 
[S II]/Ha -|-(0.083 ± 0.065) with rms = 0.081, if one 


outlier is excluded). This strong correlation shows that 
[N II]/Ha can be also useful to selecting SNRs, when 
[S II] lines are not available. Third, the SNRs in other 
galaxies also show similar correlatons, but they are all 
located below the sequence for M81 SNRs, with larger 
scatters. 

In the study of M33 SNRs, IGordon et al.l (1199811 noted 
that [O I]/Ha ratios show a strong correlation with 
[S II]/Ha ratios. In Figure lTTT bl we display the relation 
between the [O I]A6300/Ha and [S II]/Ha ratios of the 
SNRs in M81 in comparison with M31 SNRs and M33 
SNRs. M81 SNRs also show a reasonable correlation be¬ 
tween these ratios. M31 SNRs also show a similar corre¬ 
lation, if seven SNRs with high [O I]A6300/Ha ratios are 
excluded. From the linear fits we derive [O I] A6300/Ha = 
(0.616 ± 0.080) [S II]/Ha -(0.384 ± 0.075) with rms = 
0.115 for M81, and [O I]A 6300/Ha = (0.664 ± 0.065) 
[S II]/Ha -(0.301 ± 0.053) with rms = 0.117 for M33. 
The slope of this relation is much smaller than the value 
for the [N II]/Ha and [S II]/Ha relations. In addi¬ 
tion, M81 SNRs are located below the M33 and M31 se¬ 
quences, which is the opposite to the case of the [N II] ra¬ 
tios. [O I]A6300/Ha ratios can be also used for SNR 
selection, when [S II] lines are not available. 

The tight correlations between [N II], [S II], and 
[O I] and little correlation between these lines and 
[O III] lines found in this study are consistent with 
the results of the modern shock models (Dopita, private 
communication). Faster shocks increase the internal UV 
radiation field, making the recombination zone more ex¬ 
tensive. Therefore, as the shock velocity increases, the 
line ratios of [N II]/Ha, [S II]/Ha, and [O I]/Ha increase, 
while [O III]/H/3 decreases. The [O III] line emitting re¬ 
gion is not in the recombination zone, so is unaffected, 
while the H/3 line that is emitted from the recombination 
zone beco mes stronger, decreasing the [O III]/H/3 ratio. 

Figure 1121 displays the [N II]/Ha ratios against the 
[S II]A6717/[S II]A6731 ratios that are an density indi¬ 
cator, for the SNRs in M81 (this study) as well as in 
M31, M33, and other disk galaxies (Galarza et al.l 119991 : 
IGordon et al.l 119981 : ILeonidaki et al. 1201311 . M81 SNRs 
show little correlation between these two ratios, as those 
in other galaxies. This shows that [N II]/Ha ratios de¬ 
pend little on the density of the SNRs. 

3.3. Comparison with AGN Classification Diagrams 

Emission line ratio diagrams are often used for 
spectral class if icatio n of AGNs as well as SNRs. 

I Baldwin et al.l (1198111 suggested classification parame¬ 
ters for the emission-line spectra of extragalactic ob¬ 
jects, arguing that the [O III]A3727/(0 III]A5007 ra¬ 
tio is efficient to distinguish photo-ionized objects 
([O III]A3727 < [O III]A5007) and shock-heated ob¬ 
jects ([O III]A3727 > [O III]A5007) (called BPT di¬ 
agrams). However, this parameter is aff e cted much 
by extinction. Later iVeilleux fc OsterbrockI (|1987ll sug¬ 
gested other parameters that depend much less on ex¬ 
tinction: [O III]A5007/H/3 versus [N II]A 6583/Ha, 
[O III]A5007/H;3 versus [S II]AA6717,31/Ha, and 
[O III]A5007/II/3 versus [O I]A6300/Ha (called VO di¬ 
agrams). Glassification boundaries on these diagrams 
were gi ven for H II regions, Seyfer t galax ies, and LIN- 
ERs b v IShields fc Filippenkol (1199011 . Later iKewlev et all 
(|2001ll presented a boundary for the starburst galaxies 
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Fig. 12.— The [N Il]/Ha vs. [S Il]A6717/[S Il]A6731 diagram 
of M81 SNRs (circles), H II region (asterisk), and M82 SNR can¬ 
didates (crosses) in this stud y, in comparison with M31 SNRs (di¬ 
amonds. IGalarza et al.|[l999lh M33 SNRs (trianglesJGordoi^t^l 
^9^, and SNRs in six disk galaxies (plusses, ILeonidaki et al. 
2013) . Open and filled circles denote [O IIl]-strong and [O III]- 
weak SNRs, respectively. 
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and AGN based on theoretical mode ls, which was re¬ 
vised for pure star-forming galaxies in iKauffmann et ahl 
(1200,81 1. These diagrams are useful also for distinguishing 
SNRs and H II regions. 

Figure 1131 displays these diagrams for M81 SNRs 
in this study. We also plotted th e data for the SNRs 
in M31 and other d isk galaxies (jGalarza et al.l 119991 
ILeonidaki et al.l[201.8[ l for comparison. Most of the SNRs 
are located above the Kauffmann’s demarcation line in 
Figure fOI al. and all of them are outside the H II region 
boundaries, while one H II region is located inside the 
H II region boundary. Most of M81 SNRs are located in 
the LINER region, while a few are in the Seyfert region. 
This shows again that M81 SNRs are shock-ionized. 

It is noted that the two M82 objects are located 
close to the boundary of the star-forming regions in the 
[O III]A5007/H/3 versus [N II]A 6583/Ha diagram, and 
are inside the boundary in the [O III]A5007/H/3 ver¬ 
sus [S II]AA6717,31/Ha, and [O III]A5007/H/3 versus 
[O I]A6300/Ha , as some of the [O III]-weak objects are. 
In Figure [lH we plotted the same diagr ams as Figure fCT 
overlaying shock ionization models that iHo et all (120141) 
presented for 12-|-Log(0/H)=9.14, shock fraction from 0 
to 100%. and s hock velocity from 100 to 300 km s“^. 
iHo et al.l (120141) used these models to explain the shock 
features in the outflow region of star-forming galaxies 
in terms of combination of shock and photo-ionization. 
These models were given for a fixed oxygen abundance 
value, but they are still useful as a reference to explain 
the data for M82. The data of the M82 objects are con¬ 
sistent with the fast shock models with a low fraction of 
shock in the figure. This and their location in the galactic 
outflow region indicate that the two M82 objects may be 
shock condensations in the general outflow perpendicular 




Log([01] X6300/Ha) 


+ Disk SNRs (Leonidaki+l3) 
O M31 SNRs (Galarza+99) 

O M81 SNRs 
X M82 SNRs 
X M81 HII region 


Fig. 13.— Spectral classification diagrams (BPT or VO di¬ 
agrams) in the logarithmic scales of line ratios of the SNRs 
in M81 (circles, this study) in co mparison with those in 
M31 (diamo nds. IGa larza et al.| |1999l) and other disk galax¬ 
ies (plusses, ILeonidak^^T.1 12013) : (a) [O IIl]A5007/H/3 vs. 

[N Il]A6583/Ha, (b) [O IIl]A5007/H/3 vs. [S Il]AA6717,31/Ha, and 
(c) [O IIl]A5007/H/3 vs. [O l]A6300/Ha. The curved solid line and 
dashed line denote a theoretical upper limit for starburst galaxies 
(IKewlev __et__aLII2QQlJ). and a boundary for pure star-forming galax- 
ies UKauffmann et al.|[2003) . Note that all 26 targets in M81 are 
located in the SNR region and that M81 SNRs are located in two 
groups: an [O IIl]-strong group (open circles) and an [O IIl]-weak 
group (filled circles). Two M82 objects are located around the 
boundary between H II regions and SNRs, while one M81 H II re¬ 
gion is inside the H II region boundary. 



Log([NII] ^6583/Ha) Log([SII] ?.?.67I7,31/Ha) Log([OI] >.6300/Ha) 


Fig. 14.— Sam e as Fig. 13, b ut for shock ionization models 
(crosses), given in IHo et al.l II2014I) . for 12-|-Log(0/H)=9.14, shock 
fraction from 0 to 100 % with 20 % interval (from bottom to 
top), and shock velocity from 100 to 300 km s“^ (with 20 km 
s“^ interval for shock velocity 100 to 200 km s“^, and 25 km 
s“^ interval for shock velocity 200 to 300 km s“^) (from left to 
right). 


to the major axis of M82. However, the shell structures 
seen in Figure 4 indicate that they may be SNRs as well. 

3.4. Comparison with Shock-ionization Models 
iDopita et al.l (|1984D presented the shock-ionization 
models for weak shock velocity {vg < 200 km s“^). Later 
these models were extended to the case of fast shock 
mo dels covering the shock velocity of 200 to 1000 km 
s~^ (|DoDita fc Sutherlandiri99,5l 119961 : lAllen et 'mi2008[) . 
We compare d the line rati o s of M 81 SNRs with the mod¬ 
els given by iDopita et al.l ()1984D that are useful for the 
analysis of SNRs in nearby galaxies. These models were 
applied previously to the case of M31 and M33 S NRs by 
iBlair fc KirschnCT (|1985D and iSmith et al.l (Il993f) . How- 
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ever, the models of iDopita et alJ (|1984ll are based on 
old atomic physics, which is much improved nowadays. 
Emission line ratios of shock models depend not only on 
abundances of the interstellar medium, but also on the 
complexity of partially radiative shocks, degree of mix¬ 
ing with photoionization, shock velocity, the strength of 
magnetic field, and grain destruction. However, these 
factors were not included in the old models, and no mod¬ 
ern shock model grids considering all these factors are not 
yet available. Therefore the results of determination of 
abundances of the SNRs based on the old simple mod¬ 
els should be considered only as an approximate guide. 
These results can be improved, when the new modern 
shock models are available in the future. 

We plot in Figure [IHa) the [O III]/H/3 versus 
[O II]/H/3 diagr am of M81 SNRs d erived in this study 
and M31 SNRs (jGalarza et al.lHOQ^ in comparison with 
the shock-ionization models for the shock velocity of 
Vg = 106 km s“^ a nd the fixed abund ance ratio of 
O/S = 42.8 given bv iDopita et al.l (1198411 . It is noted 
that the SNRs in M81 and M31 are located roughly 
around the model grids, but their scatter in [O n]H/3 ra¬ 
tios is much larger than the range of the grids. This 
indicates again that the large scatter in [O lI]H/3 ratios 
is mainly due to the large range in shock velocity of the 
SNRs. It is noted that the large scatter in this figure can 
also be produced if the shocks are too young to be fully 
radiative, being in partially radiative state. According 
to the radiative shock models appl ied to the cas e of th e 
microquasar S26 in NGC 7793 bv IDopita et al.l (j20I2tl . 
the [O III]/H/3 ratio decreases, while the [O II]/H,0 ratio 
increases, as shocks get older. Therefore the large scatter 
in the [O III]/H/? versus [O II]/H/3 diagram must be due 
to the combination of a large range of shock velocity and 
shock age. _ 

In Figure I15l( bl we show the [O III]/H/3 versus 
[S II]A6731 /Hq! diagram of M81 SNRs (circles) and M82 
SNRs (crosses) in comparison with the same shock ion¬ 
ization models but for varying O/S abundance ratios. 
It is seen that the [O III]-strong SNRs are located 
mostly around the grid for 0/S= 43, which is similar 
to the mean value a dopted for other grids, 0/S= 42.8 by 
iDopita et al.l ( I984f) . 

Figure ma) displays the [N II]/Ha versus 
[S II]A6731/Ha diagram of M81 SNRs in comparison 
with the same shock-ionization models for various values 
of oxygen abundance and the ratio of oxygen to nitro¬ 
gen abundance (for the abundance ratio of 0/S= 42.8). 
[S II]A6731/Ha ratios are sensitive to the ratio of oxy¬ 
gen and sulfur abundance, but loses its sensitivity at the 
high abundance. [N II]/Ha ratios are a good indicator 
for nitrogen abundance. Some of MSI SNRs are located 
beyond the model grid limits, for which abundances can¬ 
not be derived. Note also that the [N II]/Ha versus 
[S II]/Ha can be a good calibrator of the N/S abundance 
ratio, since all of these lines arise in the same region of 
the shock. 

In Figure EUb) we show the [N II]/Ha versus 
[O III]/Ha diagram of M81 SNRs in comparison with 
the same shock-ionization models. In this figure, 
[O III]/Ha ratios appear to be a good indicator of oxygen 
abundance. The [O III]-strong SNRs in MSI are located 
inside the model grid, while the [O III]-weak SNRs are 
beyond the lower limit of the model grid. These [O III]- 



Fig. 15.— The [O IIl]AA4959,5007/H/3 vs. [O II]AA3727,9/Hd di- 
agram (a) and the [O IIl]AA4959,5007/H/3 vs. [S Il]A6731/Ha dia¬ 
gram (b) of M81 SNRs (circles) and HII region (cross) in compari¬ 
son with the shock ionization model grid for shock velocity Vs = 106 
km and varying ratios of o xygen, nitrogen, carb on, and sulfur 
abundances (lines) provided bv IDopita et al.l 19841') (their Figures 
7 and 9). Open and filled circles denote [O IIl]-strong and [O III]- 
weak SNRs, respectively. 
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Fig. 16.— The [N II]AA6548,83/Ha vs. [S II]A6731/Ha diagram 
(a) and the [N Il]AA6548,83/Ha vs. [O IIl]AA4959,5007/H/3 dia¬ 
gram (b) of M81 SNRs (circles) and M82 SNRs (crosses) in com- 
parison with the s hock ionization model grid for 0/S=42.8 by 
IDopita et al.l 1119841') (their Figures 8 and 10). Open and filled 
circles denote [O Illj-strong and [O IH]-weak SNRs, respectively. 
Note that [O Hij-weak SNRs are beyond the left limit in (b). 


weak SNRs might have low abundance or low shock ve- 
lodty, or both. Thus the [O III]/Ha ratio is a better 
abundance indicator t han the [S II] A 6731/ Ha , but not 
as good as claimed bv IDopita et all l|1984ll . We derived 
nitrogen and oxygen abundances of M81 SNRs using the 
model grids in Figure 1161 listing them in Table [D 

In Figure 1171 we compare the oxygen and nitrogen 
abundances derived using the two diagrams. 0/H val¬ 
ues show a large scatter and little correlations between 
two estimates. On the other hand, N/H values show a 
correlation between the two estimates so that the esti¬ 
mated values of N/H are considered to be reliable. We 
conclude that the [N II]/Ha versus [O III]/Ha diagram 
is better for abundance estimation of the SNRs than the 
[N II]/Ha versus [S II]A6731/Ha diagram. However, it is 
noted that [O III]/Ha ratios are more affected by shock 
velocities so that we need to be cautious in using the 
oxygen abundances derived from [O III]/Ha ratios. 

We calculat ed also the va l ues o f the metallicity index 
A defined by iDopita et al.l (Il984f) . Log A = Log 0/H 
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Fig. 17.— Comparison of the oxygen (left) and nitrogen 
(right) abundances of M81 SNRs derived from the [N II]/Ho; vs. 
[S Il]A6731/Ho grid and the [N II]/Ho vs. [O III]/Ho grid. 
Dashed lines represent one-to-one relations. 


+ Log N/H + Log S/H, for the fixed value of 0/S= 
42.8. These values are listed in Tabled The metallicity 
index A can be a good metallicity indicator, if reliable 
values of oxygen, nitrogen, and sulfur abundance can be 
derived from the data of SNRs. However, oxygen lines 
can be affected significantly by shock conditions, while 
nitrogen lines depend mainly on abundance. From this 
we consider that the metallicity index A derived in this 
study is not a good indicator of metallicity for SNRs. 
Therefore we do not discuss the values of A derived in 
this study, although we list them in the table for future 
studies. 

3.5. Relations between Line Ratios and Sizes of SNRs 

In Figure 1181 we plot the diameter versus galactocen- 
tric distance of M81 SNRs (this stu dy) in comparison 
with M33 SNRs (|Gordon et al.l [19^ . It is noted that 
four of six large SNRs with D > 70 kpc are [O III]- 
weak SNRs. It appears that M81 SNRs show a weak 
correlation between the size and the galactocentric dis¬ 
tance. However, a linear fit for the sample of all SNRs 
in M81 yields a slope of dD/dR = 2.97 ± 3.58, and the 
correlation coefficient is as small as 0.22, showing little 
correlation between the two parameters. Similarly M33 
SNRs show little cor relation between the two parameters 
(|Gordon et aklflO^ . It is noted that only small SNRs 
with < 60 pc are seen in the inner region at i? < 4.5 
kpc, while a large range of SNRs are seen in the outer 
region at i? > 4.5 kpc. This trend is also found for the 
SNRs in M31 and M33 (|Lee fc Leell2014al fO). This trend 
is consistent with the explanation that the pressure in the 
interstellar medium is higher closer to the galaxy center, 
and the SNRs become radiative more quickly in the inner 
region by consequence so t hat they evolve faster and die 
away (ibonita et al. l[20Toli . 

Figure 1191 displays the relations of four line ratios 
([N II]/Ha, [S II]A6717/[S II]A6731, [S II]/Ha, and 
[O I]/Ha) and sizes of M81 SNRs (this stud y) in com¬ 
parison with M33 SNRs (|Gordon et al.l[l9^ . This fig¬ 
ure shows a few notable features. First, small SNRs 
with £) < 70 pc in M81 show a larger scatter in the 
ratios of these lines. Second, large SNRs with D > 70 
pc are mostly (four out of six) [O Illj-weak SNRs. The 
[O III]-weak SNRs show lower values with a much smaller 
scatter in the ratios of these lines. However, two large 
SNRs (MF 06, MF 17) have high values. Considering 
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Fig. 18.— The size vs. galactocentric distance of M81 SNRs 
(circles) in this stu dy in comparison with M33 SNRs (triangles, 
IGordon et al.lfT998l) . Open and filled circles denote [O Hl]-strong 
and [O Hl]-weak SNRs, respectively. Dashed line represent a linear 
fit for the sample of all SNRs in M81, with a slope of 2.97di3.58. It 
appears to show a weak correlation between the size and galacto¬ 
centric distance, but the correlation coefficient is as small as 0.22, 
showing little correlation. 
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Fig. 19.— Emission line ratios vs. size of M81 SN Rs (circles) in 
this study in comparison with M33 SNRs (triangles. IGordon et all 
[19^ : (a) [N Il]/Ha , (b) [S Il]A6717/[S Il]A6731, (c) [S Il]/Ha, 
and (d) [O l]A6300/Ha. Open and filled circles denote [O IIl]- 
strong and [O Hl]-weak SNRs, respectively. 
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Fig. 20. — The [S Il]A6717/[S Il]A6731 ratio vs. galactocentric 
distance of M81 SN Rs (cir cles) in this study with M33 SNRs (tri¬ 
angles, [Gordo^^^3]ES3)- Open and filled circles denote [O III]- 
strong and [O IIl]-weak SNRs, respectively. 

that the [O III]-weak SNRs have fast shock velocity and 
they are relatively large, we conclude that they may be 
dynamically old. Third, small and large SNRs in M33 
show similar behaviors. However, M81 SNRs have much 
higher [N II]/Ha ratios than M33 SNRs. Fourth, the 
[S II]A6717/[S II]A6731 ratios of small SNRs in M81 (as 
well as in M33) are on average larger than those of the 
large SNRs, showing that larger SNRs have in general 
lower densities than smaller SNRs, as expected. 

3.6. Radial Variation of Sizes, Line Ratios, and 
Abundances of SNRs 

Figure [20] plots the [S II]A6717/[S II]A6731 ratio ver¬ 
sus ga lactocentric distanc e of M81 SNRs as well M33 
SNRs ([Gordon et al.lll998ll . What is the most impressive 
about this figure is that M81 SNRs show clearly an upper 
envelope at [S II]A6717/[S II]A6731 ss 1.45. In compari¬ 
son, M33 SNRs show a larger scatter, and some of them 
have [S II]A6717/[S II]A6731 values above the envelope. 
The latter might be due to the larger errors in the mea¬ 
surements of M33 SNRs ([Gordon et al.lIlQQSll . The value 
for the upper envelope of M81 SNRs is close to a theo¬ 
retical lower limit, 1.43, ([Blair fc Kirschn^ll985[l . This 
shows that the measurements of [S II]A6717/[S II]A6731 
ratios of M81 SNRs in this study are solid. All M81 SNRs 
have [S II]A6717/[S II]A6731 > 1.15. This value corre¬ 
sponds to the density of 312 cm“^, according to the con¬ 
version of [S II]A6717/[S II]A6731 ratio into densities as a 
functi on of temperature in Figure 7 of iBlair fc Kirschned 
([1985fl . This shows that they are located in the low den¬ 
sity region. M81 SNRs show little radial gradient in this 
line ratio. 

We investigated any radial gradient of emission line ra¬ 
tios of M81 SNRs. In Figure 1211 we displayed four line 
ratios ([O III]A5007/H/3, [O III]A5007/[0 II], [N II]/Ha, 
and [S II]/Ha ratios) versus galactocentric distance of 




Fig. 21.— Several line ratios vs. galactocentric di stance for M81 
SNRs in comparison with M31 S NRs (diamonds. IGalarza et ahl 
1199911 and M33 SNRs (triangles, I Gordon et al.l 1199811 : (a) Log 

[O IIlJ/Ha, (b) Log [O IIl]/[0 II], (c) Log [N IlJ/Ho , and (d) 
Log [S Il]/Ha. Open and filled circles denote [O IIl]-strong and 
[O IIl]-weak SNRs, respectively. The solid line and dashed line in 
(c) represent linear fits for all M81 SNRs and for M81 SNRs ex¬ 
cluding one outlier at about 12 kpc, with slopes, —0.013 ih 0.007 
and —0.018 zb 0.008 dex kpc“^, respectively. The do tted-line is a 
linear fit for M31 SNRs given bv IGalarza et ahl 1)199911 with a slope 
—0.04 dex kpc“^, and the dot-dashed line a linear for M33 SNRs, 
with a slope —0.05 dex kpc“^. The solid-line in (d) represents a 
linear fit for all M81 SNRs with a slope, —0.015 ih 0.008 dex kpc“. 


M81 SNRs. We plotted also the data for M31 and 
M33 ([Galarza et ahl 119991 : iGordon et ahl I1998D . Several 
features are noted in this figure. First, M81 SNRs 
show little radial gradients in [O III]A5007/H/3 and 
[O III]A5007/[0 II] ratios. Second, M81 SNRs show 
clearly a radial gradient in [N II]/Ha ratios. Linear 
fitting for 26 SNRs yields Log [N II]/Ha = (—0.013 ± 
0.007)R -b (0.085 ± 0.048) with rms = 0.09 (^0.018 ± 
0.008)i? + (0.113 ± 0.055) with rms = 0.09 if two out¬ 
liers are excluded). Even if we select only the [O III]- 
strong SNRs, we obtain similar results: Log [N II]/Ha = 
(-0.015 ± 0.006)R -b (0.132 ± 0.040) with rms = 0.07. 
Similarly we derive a linear fit for 72 SNRs in M33: Log 
[N II]/Ha = ( -0.050 ± 0.006)R - ( 0.238 ± 0.026) with 
rms = 0.124. IGalarza et al.l (|1999H presented a value 
of the radial gradient for M31 SNRs, —0.04 ± 0.01 dex 
kpc“^. Thus the slope of M81 SNRs is two to three times 
flatter than those of M31 SNRs and M33 SNRs. Third, 
the [S II]/Ha ratios of M81 SNRs also show a radial gra¬ 
dient. Linear fitting yields Log [S II]/Ha = (—0.016 ± 
0.008)i?-b(0.104±0.057) with rms = 0.09, showing a sim¬ 
ilar slope to that of the [N II]/Ha ratios. Even if only the 
[O III]-strong SNRs are used, we obtain similar results: 
Log [S II]/Ha = (-0.017±0.008)R-b(0.111 ±0.051) with 
rms = 0.08. This is consistent with the tight correlation 
between these two line ratios, as described before. 

In Figure 1221 we plotted Log(N/H) and Log(0/H) 
versus galactocentric distance of M81 SNRs derived 
from the comparison with the [N II]/Ha versus 
[O III]A5007/H/3 grid of the shock ionization models 
in the previous section. We also plotted the data 
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Fig. 22.— Nitrogen (a) and oxygen (b) abundance vs. galac- 
tocentric distance of M81 SNRs (circles) in this stu dy, in com¬ 
parison with M33 SNRs fnlusses. ISmith et al.lfT993l’) . Note that 
the abundances of M81 SNRs were derived from the [N II]/Ha vs. 
[O III]/Ha grid in this study, while thos e of M33 SNRs wer e from 
the [N H]/Ha vs. [S H]A6731/Ha grid bv ISmith et al.l 1119931) . 


for M33 SNRs derived from the [N II] /Ha versus 
[S II]A 6 731/Hq! grid of the same models in iSmith et al.l 
(|1993ll . ISmith et al.l (|1993D covered only the red wave¬ 
length in their spectra so that they had no data for 
the [O III] lines and could not use the [N II]/Ha ver¬ 
sus [O III]/Ha grid for their analysis. Fortunately M33 
SNRs have low abundance so that they could estimate 
the abundance of M33 SNRs using the [N II]/Ha ver¬ 
sus [S II]A6731/Ha grid. M81 SNRs show clearly a ra¬ 
dial gradient in Log (N/H), and we derived from lin¬ 
ear fitting for 21 SNRs at 2 < i? < 14 kpc, Log 
(N/H) +12 = (-0.023 ± 0.009)R -f (8.154 ± 0.069) with 
rms = 0.184 (= (-0.020±0.007)R-f (8.101 ±0.049) with 
rms = 0.125, if two outliers are excluded). This slope 
is much flatter than that of M3 3 SNRs (0.5 < R < 4.5 
kpc), —0.14 ± 0.04 dex kpc“^ (jSmith et al.lllQQ.l l. On 
the other ha nd, M81 SNRs show little radial gradient in 
Log (0/H). ISmith et al.l 1)19931) presented a weak oxy¬ 
gen gradient for M33 SNRs, —0.09 ±0.05 dex kpc“^, but 
with a large scatter. These results for the abundances 
are consistent with the results for the line ratios for M81 
SNRs. 

It is known that the pre-SN WR stars can enrich their 
surroundings with CN processed gas from mass-loss, as 
seen i n the study of WR ring nebulae such as NGC 
6888 (jMesa-Delgado et 1111 1201^ . This would raise the 
N abundance, while keeping S and O abundances un¬ 
changed. However, M81 SNRs show radial gradients in 
both [N ll]/Ha and [S II]/Ha ratios as in Figure 19. 
Therefore the contribution of this effect is considered to 
be not significant for the nitrogen gradient of M81 SNRs. 
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Fig. 23.— 12-l-log(0/H) vs. galactocentric distance for the 
H H regions in M81 and two M82 objects (red filled circles) 
in this study i n comp arison with the M81 H H regions from 
IPatterson et al.1 II2Q12I) (blue open circles). The outermost two 
objects belong to M82. Oxygen abundances were derive d us- 
ing four different calibration s: (a) lPilyugin_^£_Thuan| ([200^ ), (b) 
IKobulnickv &: KewlevI (120041). (c) iKewle^J^Comta 1120021'). and 
(dju res^ 11III2007I). 


We derived the physical parameters of one H II region 
in M81 using the strong-line analysis method, following 
the descriptio n given in the manual fo r the Nebular pack¬ 
age in IRAF ()Shaw fc Dufourl [199511 . We also applied 
this procedure to two M82 objects, in case that they 
are H 11 regions, rather than SNRs. We used four em¬ 
pirical and theoretical calibrations for deriving oxygen 
abundances from emission line fluxes of these objects , 
as done for M81 H 11 re gions in IPatterson et ahl (|2012f l 
and M31 HH re gions in Sanders_etaL (120121) : empiri- 
cal cal ibrations in iPilvugin fc: Thuanl (1200511 andlBresolin 
(2007), an d theoretical calibrations iniKewlev fc Dooita 
( 2002 ) and IKobulnickv fc KewlevI (|2004D . Table l5] lists a 
summary of the oxygen abundance values derived in this 
study. _ 

Figure 1231 shows the oxygen abundances of one 
M81 H 11 region we derived versus the galactocen¬ 
tric distance, in co mparison with those of other M81 
H II regions given by IPatterson et al.1 (|2012[1 . The value 
of our H II region i s consistent with the results of 
IPatterson et al.1 (120121) for the same galactocentric dis¬ 
tance. On the other hand, two objects in M82 show 
the oxygen abundance s higher than the mean value of 
the M81 H II regions. ISmith et al.1 (|2006f ) derived from 
HST/STIS spectroscopy a value of the oxygen abun¬ 
dance for one compact H II region hosting a star clus- 
te r M82 A-1: Log(0/H)±12=9.23, 8.76, and 8.77, based 
on iKewlev fc Dopital (120021) . IPettini &: Pagell (|2004l) . and 
adoption of a value of Tg = 10^ K, respectively. These 
results are close to or higher than the solar value, 
Log(0/H)=8.7 to 8.9. Thus the oxygen abundances of 
the two M82 objects derived in this study are similar to 
the value for M82 A-1. 


4. DISCUSSION 

4.1. Comparison of Abundance Gradients of SNRs, 
HII Regions, and PNe in M81 


3.7. Results for HII Regions 
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We compared the radial gradients of nitrogen and oxy¬ 
gen of SNRs derived in this study with those based 
on H II regions and PNe in M81 in the literature 
(jPatterson et, al.ll2 012tlStanghellini et al.ll2ni^ . as listed 
in Table [S] IPatterson et all (|2012ll analyzed a sample 
of H II regions at the large range of 5.7 < R < 32 kpc 
in M81, using the method of strong-line oxygen abun¬ 
dance analysis. They derived a value for the radial oxy¬ 
gen gradient, —0.013 to —0.020 dex kpc“^, much flatter 
than the previous results based on the samples of H II re¬ 
gions in the inner region of M81 (IGarnett fc Shieldill987t 
IStanghellini et al.]l2010ll. _ 

On the other hand. Istanghellini et al.l (|2014f) derived 
abundances of H II regions and PNe in the inner re¬ 
gion {R < 13 kpc) of M81, using the method of weak- 
line abundance analysis. They suggested that the ra¬ 
dial oxygen gradient of H II regions in the inner re- 
gion (R < 10 k pc) o f M81 is much steeper than the 
IPatterson et al.l (1201211 ’s value. Then they pointed out 
that there may be a break at i? ~ 15 kpc, noting 
that two H II regions at 16 kpc and 32 kpc given by 
IPatterson et al.l (|2012ll showed a constant value at lower 
oxygen abundance, log (0/H)-|-12 « 8.1 (see their Fig. 
8). They found also that the radial oxygen gradient of 
the H II regions at 4 < i? < 10 kpc (—0.088 ± 0.013 
dex kpc“^) shows an about twice steeper slope than the 
radial metallicity gradient of PNe at 3 < i? < 13 kpc 
(—0.044 ± 0.007 dex kpc“^), suggesting that this is due 
to evolu tion effect of meta l licity gradient (see their Fig¬ 
ure 5). Istanghellini et aH (|2014ll found that the radial 
nitrogen gradient of the H II regions (—0.067 ± 0.013 dex 
kpc“^) shows an about 1.5 times steeper slope than the 
radial metallicity gradient of PNe (—0.049 ± 0.007 dex 
kpc“^) (see their Figure 7). They also pointed out that 
these gradients in M81 are steeper than those in other 
spiral galaxies (MWG, M33, and NGG 300), and that 
this may be due to its location in the galaxy group (see 
th eir Figure 9). The di screp ancy between the results 
of IPatterson et al.l (|2012ll and IStanghellini et al.1 (|2014f) 
may be due to the differences in the method used and 
the radial coverage of the samples in the two studies, 
and further studies of H II regions in M81 are needed to 
clarify this problem. 

The radial gradient of nitrogen abundance of the SNRs 
at 2 < i? < 14 kpc in M81 derived in this study, —0.023± 
0.009 dex kpc“^, is a few times flatter than the values for 
the nitrogen abundance gradient of the H II regions and 
PNe g iven by Istanghellini et ahl (120141) . IPatterson et al.l 
(I2012fl did not present any data for the nitrogen radial 
gradient of the H II regi ons. However, the flatt er oxygen 
radial gradient given bv IPatterson et all (|2012ll indicates 
a similarly flatter nitro gen radial gradient , much flatter 
than the values given bv Istanghellini et al.l (120141) . Then 
little radial gradient of oxygen abundance of M81 SNRs 
(—0.010 ±0.015) found in this study may be close r to the 
flatter gradient given by IPatterson et al.l (120121) r ather 
than t o the steeper gradient given bv istanghellini et all 

(IMl . 


Why SNRs show flatter radial gradients in nitrogen 
and oxygen than H II regions and PNe in M81 is not 
clear, requiring further studies. SNRs located in the spi¬ 
ral arms as those in this study mostly must have come 
from core-collapse supernovae in star-forming regions. 


Then they are as young as H II regions. From this we 
expect that SNRs and H II regions follow similar radial 
abundance gradients. However, optical emission lines 
from SNRs depend both on abundance and shock condi¬ 
tions, while those from H II regions depend on abundance 
only. Oxygen lines in SNRs are more affected by shock 
conditions than nitrogen lines. Then it is expected that 
radial abundance gradients derived from optical spectra 
of SNRs and H II regions will be similar for nitrogen, 
but different for oxygen. However, the radial gradient 
for nitrogen of SNRs in this study is significantly flat¬ 
ter th an that of H II regions given by Istanghellini et aH 
(|2014f) . The cause for this discrepancy may be due to 1) 
the uncertainty in deriving abundance from the compar¬ 
ison of emission lines and shock ionization models, or 2) 
any evolution effect related with nitrogen. 


4.2. Comparison of Abundance Gradients of SNRs in 
Nearby Galaxies 

In Table in we list previous estimate s of the abundance 
gradi ents b ased on SNRs in M31 (iBlair fc Kirschii^ 
1 19851 ). M33 (|Blair fc Kirschne^ll985HSmith et iahl 1993D. 
the Milky Way Galaxy (MWG) (|Blair fc KirschnCT 
119851) . as well as in M81 in this study. M31 SNRs and 
M81 SNRs are similar in that they show a radial gradient 
in nitrogen abundance, but little in oxygen abundance, 
although the nitrogen gradient of M31 SNRs (-0.04) is 
twice steeper than that of M81. On the other hand, 
M33 SNRs show radial gradients in both nitrogen and 
oxygen, which are steeper than those for M81 and M31. 
The MWG SNRs show a weak radial gradient of nitrogen, 
similar to M81 SNRs. There is no information available 
fo r the oxygen radial gradi ent of the SNRs in the MWG. 

IBlair fc Kirschn^ (| 198 51 ) noted that the SNRs in M31 
and M33 show little radial gradients and a large scat¬ 
ter in oxygen abundances, while they show clearly radial 
gradients similar to H II regions but a much larger mean 
values in nitrogen abundance. They pointed out that the 
causes for the oxygen abundance are two: the possible 
confusion with low shock velocity SNRs, and the con¬ 
tamination due to nearby photoionized regions. In the 
case of M81 SNRs in this study, we could separate the 
SNRs with low and high shock velocities according to the 
[O IIIJA5007/H/3 ratios, [O III]-weak and [O III]-strong 
groups. Even if we select only the [O III]-strong SNRs, 
they show little radial oxygen gradients. As described in 
the previous section, the value of the [O III]A5007/H/3 ra¬ 
tio changes rapidly from -0.5 at Vg « 80 km s“^ to 
±0.5 at Us « 100 km s“^, and increases slowly to 
ss 0.65 at Vg ^ 170 km s“^, according to the model of 
iDopita et al.l (|1984D . On the other hand, the values of 
the [O IIIJA5007/H/3 ratio for the [O III]-strong SNRs in 
M81 derived in this study are 0.0 to 0.7. Therefore little 
radial gradients and a large scatter in oxygen abundances 
in M81 SNRs are mainly due to the varying shock veloc¬ 
ity even among the [O III]-strong SNRs. The spectra of 
M81 SNRs were obtained with fibers with 1'.'5 diameter 
so that it is possible that they might have been contami¬ 
nated by nearby photoionized regions. However, most of 
M81 SNRs in this study are larger than the fiber sizes so 
that the effect of contamination due to nearby photoion¬ 
ized regions is not considered to be significant. 


4.3. Optical SNRs and X-ray Sources in M81 
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Fig. 24.— X-ray hardness r atio diagram (X-ray classification di¬ 
agram, [Pre3wic^^^23[2SS3) of the optical SNRs (filled starlets) 
matched with the catalog of M81 X-ray sources (open diamonds, 
ISell et al.lf2QTI^ . S, M, and H represent, respectively, the counts 
in soft (0.5-1 keV), medium (1-2 keV), and hard (2—8 keV) bands, 
and T for the total counts (0.5-8 keV). Dotted lines denote the 
boundary for net counts. Note that the positions of the optical 
SNRs are consistent with the region for thermal SNRs (large cir¬ 
cle). The X-ray sources with harder colors are mostly HMXBs and 
LMXBs. 

iPanmiti et alJ (j2007fl found 97 X-ray sources in M81 
through the Chandra X-ray Observatory (CXO) obser¬ 
vations. They searched for X-ray counte rparts of the op¬ 
tical and radio SNR candidates in M81 (|Matonick et al.l 
1997li using th is catalog, but finding none. Later 
Sell et al] (|2011ll provided a catalog of 265 X-ray sources 
in M81 detected in a large number of CXO fields, which is 
also included in iLfi] (prill . We crosschecked our sample 
of M81 S NRs with t he ca talog of X-ray sources in M81 
given by ISell et all l) 201 111 , finding five SNRs matched 
with X-ray sources: MF10=Sell259, MF17=Selll93, 
MF19=Selll95=Liul578, MF25=Selll72, and 
L5=Sell50=Liul653. The images of these SNRs in 
Figure [3] show that they are compact and bright, 
indicating that they are relatively young. They are 
all located in the inner spiral arm region, as shown in 
Figure 1. They all have [O > 1, belonging to 

the [O III]-strong group, which show that they have 
relatively high shock velocity. 

In Figure [24] we plot these SNRs in the X-ra y hardness 
ratio diagram of M81 X-ray sources given by ISell et all 
(Unni). S, M, and H represent, respectively, the counts 
in soft (0.5-1 keV), medium (1-2 keV), and hard (2- 
8 keV) bands, and T for the total counts (0.5-8 keV). 
(M — S)/T and (iJ — M)/T denote soft X-ray color and 
hard X-ray color, respectively. This diagram is useful 
for classification of the X-ray sources in nearby galax¬ 
ies to distinguish Low Mass X-ray Binaries (LMXBs), 
High Mass X-ray Binari es (HMXBs), SNRs, absorbed 
AGN sources and stars ([Prestwich et al.l [200311 . X-ray 
sources in the high density region with hard colors are 
mostly HMXBs and LMXBs that are composed of either 
pulsar binaries or blackho le binaries and low rn agnetic 
field neutron star binaries (jPrestwich et al.ll2003ll . Ther¬ 
mal X-ray radiation is emitted from the hot plasma in¬ 
side the SNRs, while non-thermal X-ray rad iation come s 
from the shocked region in the shell regions (|Vin innii). 


The contribution of thermal line emission is strong in the 
medium band so that the thermal SNRs are expected to 
be located in the soft color region as shown in this fig¬ 
ure. Some of the X-ray sources located in the SNR region 
show variability. They a re not SN Rs. because SNRs do 
not show any variability (ILiull201lil . 

All five SNRs in M81 have soft colors and they are lo¬ 
cated in the thermal SNR region in the figure. This shows 
that these SNRs are all thermal SNRs. Thus these SNRs 
are genuine SNRs in both optical and X-ray properties. 
One of them has (M — S)/T ss —1.0, which is similar to 
the colors of the supersoft sources. The possible origins 
of the supersoft sources are either SNRs or accretion- 
powered sources that are variabl e and may be associate d 
with young stellar population ([Prestwich et al.l l2003f) . 
The presence of one optical SNR with supersoft colors 
in M81 provide support for the SNR origin of the super- 
soft sources. 

5. SUMMARY AND CONCLUSIONS 

We obtained optical spectra of 26 SNRs and one 
H II region in M81 and two SNRs in M82, covering a 
wide range of wavelength. Analyzing the emission line 
fluxes of these objects in comparison with shock ioniza¬ 
tion models for SNRs, we investigated three aspects of 
the relations: between line ratios, between line ratios and 
sizes, and between line ratios as well as abundance and 
galactocentric distances. We also derived oxygen abun¬ 
dance of one H II region in M81 and two objects in M82 
through strong-line analysis. Main findings are as fol¬ 
lows. 

1. Twenty six out of the MSI SNR candidates turn 
out to be genuine SNRs, showing [S II]/Ha ratios 
larger 0.5. Two objects in M82 are shock conden¬ 
sations in the general outflow or SNRs. MSI SNRs 
are mostly located along the spiral arms. Their 
kinematics follows well the rotation of the M81 
disk, showing that all of them belong to the disk 
population. 

2. The distribution of [N II]/Ha ratios of MSI SNRs 
is bimodal. The high component is mainly due to 
the presence of the SNRs at the inner region (i? < 5 
kpc). 

3. M81 SNRs are divided into two groups in the shock- 
ionization regions of the spectral line ratio dia¬ 
grams: an [O Illj-weak group ([O III]/H/3 < I) 
and an [O III]-strong group ([O III]/H/3 > I). The 
[O Illj-weak SNRs may have fast shock velocity 
and a low fraction of shock, and they are relatively 
large. This implies that they may be dynamically 
old. 

4. [N II]/Ha and [S II]/Ha ratios of MSI SNRs show 
a strong correlation. These ratios can be a good 
calibrator of N/S abundance ratio, because these 
lines are coming from the same region of the shock. 
[O I]/Ha and [S II]/Ha ratios also show a similar 
correlation, but with a larger scatter. 

5. [S II]A6717/[S II]A6731 ratios of the SNRs are 
smaller than 1.45, consistent with the theoretical 
density lower limit. Most of them are larger than 
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1.3, indicating that the SNRs are in the low density 
region. 

6. [N [S and [S II]A6717/[S II]A6731 

ratios of the SNRs show little dependence on the 
SNR size. 

7. [N II]/Ha and [S II]/Ha ratios of the SNRs show 

a clear radial gradient, dLog ([N II]/HQ;)/dLog 
R = —0.018 ± 0.008 dex kpc“^ and dLog 

([S II]/Ha)/dLog R = —0.016 ± 0.008 dex kpc“^ 
for 2 < i? < 14 kpc. 

8. We estimated the nitrogen and oxygen abundance 
of the SNRs from the comparison of their line ratios 
with shock-ionization models. Note that the results 
of determination of abundances of the SNRs based 
on the old simple models should be considered only 
as an approximate guide. We found a value for 
the nitrogen radial gradient, dLog(N/H)/dLogR = 
—0.023 ± 0.009 dex kpc“^, and little gradient for 
oxygen. Little gradient of oxygen of SNRs is mainly 
due to the fact that oxygen lines are significantly 
affected by the shock velocity and M81 SNRs have 
a large range of shock velocity. 

9. The nitrogen abundance of the SNRs shows a few 
times flatter gradient than those of the H II regions 
and PNe in M81. The difference in the radial gradi¬ 
ents between SNRs, PNe, and H II regions in M81 
remain to be explained. 

10. We found X-ray SNRs in M81 for the first time. 
Five of the SNRs in M81 that were confirmed from 
optical spectra are matched with X-ray sources de¬ 
tected in the CXO observations. They look mostly 
compact and bright in the [S II] images. They are 
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located in the inner spiral arm region at i? < 8.3 
kpc. Their X-ray hardness colors are consistent 
with thermal SNRs. The X-ray hardness colors for 
one of them are similar to those of the supersoft 
X-ray sources, supporting the SNR hypothesis for 
the origin of the supersoft sources. 


Through this study based on high quality spectra we 
found several interesting features of SNRs and SNR can¬ 
didates in M81 and M82, providing strong clues to un¬ 
derstand the statistical properties of SNRs. Emission 
line ratios of shock models for SNRs depend not only on 
abundances of the interstellar medium, but also on the 
complexity of partially radiative shocks, degree of mix¬ 
ing with photoionization, shock velocity, the strength of 
magnetic field, and grain destruction. However, these 
factors were not included in the old models used in this 
study, and no modern shock model grids considering all 
these factors are not yet available. Therefore the results 
of determination of abundances of the SNRs based on 
the old simple models should be considered only as an 
approximate guide. These results can be improved, when 
the new modern shock models are available in the future. 
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TABLE 1 

Emission-line Strengths oe Emission-line Objects in M81 and M82 


ID 

c(H/3) 


[O II]3727,9‘’ 

[O III]4959‘> 

[O III]A5007‘= 

[O I] 6300'^ 

[N II]6548‘> 


[N iijesssi’ 

[S II] 6717'^ 

[S II]673l‘> 

LOl 

0.42 ± 0.04 

5.9 ±0.6 

380.2 ±20.3 

17.9 ±2.6 

54.3 ±3.0 

17.8 ± 3.9 

58.7 ± 1.9 

321.6 ± 10.2 

175.1 ±5.8 

138.7 ± 5.1 

100.2 ±3.9 

L02 

0.41 ± 0.06 

3.8 ±0.5 

390.6 ±26.1 

9.4 ±3.7 

27.4 ±3.4 

20.9 ±7.1 

61.6 ±2.6 

321.1 ± 13.1 

184.1 ±7.9 

121.2 ±5.9 

90.6 ±4.8 

LOS 

0.59 ±0.07 

0.3 ±0.1 

368.0 ±31.2 

50.6 ±4.7 

131.6 ±5.0 

14.9 ±5.0 

68.7 ±3.3 

330.5 ± 15.7 

204.8 ± 10.0 

124.2 ±6.8 

93.7 ± 5.5 

L04 

0.46 ± 0.04 

3.6 ±0.3 

422.6 ± 17.0 

25.7 ± 2.6 

58.6 ±2.5 

34.8 ±4.4 

63.5 ± 1.7 

323.4 ± 8.4 

189.5 ±5.1 

142.6 ±4.3 

98.8 ±3.2 

LOS* 

0.47 ± 0.03 

5.8 ±0.4 

674.6 ± 19.1 

106.8 ±2.3 

321.0 ±3.1 

57.4 ±2.5 

98.4 ±2.1 

324.2 ±6.9 

293.7 ±6.2 

175.0 ±4.2 

151.1 ±3.7 

L06 

0.21 ±0.05 

12.9 ± 1.7 

393.1 ±25.0 

91.0 ±4.4 

262.4 ±5.4 

60.0 ±8.7 

54.7 ±2.4 

310.4 ± 12.6 

163.7 ±7.1 

102.5 ±5.3 

73.0 ± 4.2 

MFOl 

0.45 ± 0.07 

1.5 ±0.2 

492.9 ±35.6 

110.7 ±5.4 

336.4 ± 6.7 

34.6 ± 10.0 

72.0 ±3.6 

323.0 ±15.7 

215.0 ± 10.8 

153.7 ±8.5 

108.0 ±6.4 

MF04 

0.59 ±0.11 

0.5 ±0.1 

511.9 ±58.1 

55.9 ±8.3 

149.0 ±8.6 

38.6 ± 16.4 

74.7 ±5.7 

330.6 ±25.1 

222.7 ± 17.4 

172.4 ± 14.6 

137.4 ± 12.1 

MF05 

0.38 ± 0.08 

0.9 ±0.2 

601.3 ±44.9 

44.4 ± 5.9 

134.1 ±6.3 

69.6 ±9.7 

62.5 ±3.6 

319.4 ±17.6 

186.6 ± 10.8 

141.8 ±9.1 

98.7 ± 7.0 

MF07 

0.62 ±0.11 

0.5 ±0.1 

1103.0 ± 107.4 

64.8 ±9.1 

197.1 ±9.4 

149.1 ± 20.1 

88.4 ±7.2 

332.4 ±26.9 

263.4 ±21.6 

220.3 ± 19.4 

160.2 ± 14.6 

MF08 

0.74 ±0.17 

0.0 ±0.0 

1549.1 ±221.0 

132.8 ±17.0 

282.9 ± 15.7 

40.7 ± 11.1 

121.1 ± 14.7 

338.7 ±42.0 

359.4 ±44.5 

252.0 ±33.6 

180.9 ±25.1 

MFIO* 

0.67 ±0.07 

0.9 ±0.1 

877.4 ±52.4 

75.0 ±5.0 

233.8 ±5.9 

107.8 ±9.3 

95.3 ±4.5 

335.2 ± 15.9 

283.9 ± 13.7 

202.0 ± 10.6 

150.7 ±8.3 

MFll 

1.01 ±0.16 

0.1 ±0.0 

1777.0 ± 224.5 

82.2 ± 13.8 

196.0 ± 12.2 

136.1 ± 19.4 

116.9 ±12.7 

354.0 ±38.8 

346.7 ± 38.3 

261.4 ± 30.8 

182.9 ±22.1 

MF12 

0.65 ± 0.08 

0.5 ±0.1 

1214.3 ±85.2 

170.5 ± 7.3 

483.8 ±9.2 

38.4 ±8.8 

95.8 ±5.6 

333.9 ± 19.5 

285.4 ± 16.9 

198.1 ± 12.8 

139.3 ±9.5 

MF16 

0.74 ±0.11 

0.6 ±0.2 

1043.2 ±96.8 

84.2 ±9.0 

197.2 ±9.0 

47.9 ± 13.0 

76.5 ±5.8 

338.7 ±25.2 

227.8 ± 17.4 

188.6 ± 15.5 

131.0 ± 11.3 

MF17* 

0.50 ±0.04 

2.8 ±0.3 

709.6 ±28.5 

58.1 ±3.1 

167.3 ±3.6 

92.3 ±5.3 

104.7 ±3.2 

325.4 ± 10.0 

312.3 ±9.6 

234.3 ± 7.8 

172.4 ±6.0 

MF19* 

0.54 ± 0.03 

3.7 ±0.3 

655.1 ±21.6 

75.4 ± 2.5 

215.5 ±3.1 

70.4 ±4.2 

96.2 ±2.3 

327.8 ± 7.9 

286.8 ±6.9 

173.0 ±4.7 

148.9 ±4.1 

MF21 

0.34 ± 0.05 

7.6 ±0.9 

583.4 ±29.1 

37.3 ±3.2 

135.8 ±4.1 

94.3 ±5.8 

82.6 ±3.0 

317.3 ±11.4 

246.8 ±9.0 

179.2 ±7.2 

137.3 ±5.8 

MF22 

0.82 ±0.13 

0.2 ±0.0 

664.7 ±87.8 

7.4 ±6.5 

61.8 ±9.0 

113.9 ± 16.3 

78.3 ±7.0 

343.5 ± 30.6 

232.8 ±21.3 

188.9 ± 18.5 

137.6 ± 14.0 

MF25* 

0.69 ± 0.06 

0.7±0.1 

904.7 ±48.8 

91.2 ±4.7 

260.9 ±5.4 

75.1 ±7.2 

99.8 ±4.2 

335.9 ± 14.2 

297.2 ± 12.6 

194.0 ±9.0 

155.4 ±7.5 

MF26 

0.88 ±0.16 

0.0 ±0.0 

917.8 ± 133.2 

58.1 ± 15.1 

109.4 ± 13.3 

68.4 ±20.1 

102.6 ±11.6 

346.9 ± 39.2 

304.7 ± 34.8 

233.4 ± 28.5 

169.5 ±21.5 

MF27 

0.58 ±0.07 

0.4 ±0.1 

712.4 ±49.4 

101.7 ±6.0 

303.5 ± 7.2 

35.4 ±6.5 

83.5 ±4.4 

329.8 ±17.3 

248.8 ± 13.3 

164.6 ± 9.8 

124.4 ±7.8 

MF29 

0.57 ±0.06 

1.5 ±0.2 

363.8 ±28.6 

19.0 ±5.6 

27.6 ±4.3 

34.6 ±7.2 

53.8 ±2.5 

329.5 ± 14.6 

160.5 ± 7.6 

132.8 ±6.9 

91.0 ±5.1 

MF32 

0.53 ±0.10 

0.7 ±0.2 

1022.7 ±92.6 

56.5 ±8.4 

183.0 ±9.2 

150.1 ± 15.7 

93.2 ±7.0 

327.3 ±24.7 

278.0 ± 21.2 

208.7 ± 17.3 

143.7 ± 12.6 

MF33 

0.75 ±0.13 

0.6 ±0.2 

1866.7 ± 189.4 

190.2 ± 12.0 

419.5 ± 12.3 

73.9 ± 14.1 

112.4 ±9.9 

339.4 ± 30.3 

334.4 ± 30.0 

237.2 ±22.9 

182.0 ± 18.0 

MF35 

0.71 ±0.05 

3.4 ±0.4 

484.3 ±26.0 

16.6 ±3.2 

41.0 ±3.1 

108.5 ±6.1 

65.6 ±2.2 

337.1 ± 11.0 

195.4 ±6.6 

210.0 ±7.5 

156.3 ±5.8 

M82-L1 

U.bS ± U.Uvi 

1.8 ±0.1 

208.8 ±9.7 

2U.U zb l.b 

53.6 ± 1.6 

14.3 ± 1.8 

47.3 ± 1.1 

313.3 ± 8.1 

140.7 ±3.2 

59.9 ± 1.7 

46.3 ± 1.5 

M82-L2 

0.57 ±0.14 

5.8 ± 1.9 

171.8 ± 19.4 

18.7 ± 1.6 

55.1 ± 1.7 

16.9 ±2.5 

47.1 ±4.6 

307.7 ±47.1 

140.4 ±13.8 

51.1 ±5.4 

41.9 ±4.5 

M8imi-Ll 

0.46 ± 0.05 

iy.4 ± 2.2 

± 2i.U 

b5. Y ± b.b 

165.2 ± 3.8 

2.t) ± i.y 

19.7 ±0.9 

bui.y ± iu.u 

58.8 ± 2.8 

24.4 ± 2.1 

16.4 ± 1.9 


X-ray detected sources JSell et al.ll201ll) . 


^ flux in units of 10”’^^erg cm^^ s”'^ A 
^ Ail other fluxes are normalized by H/3 — 100. 
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TABLE 2 

Line Ratios and Abundances of Emission-line Objects in M81 and M82 


ID 

Size 

[pc] 

R 

[kpc] 

V 

km 

R.A. 

(J2000) 

Dec. 

(J2000) 

[OIII]A5007 

H/3 

[0116300 

Ha 

[NII16583 

Ha 

[SII16717,31 

Ha 

[SII]6717 

131116731 

log (N/H)*’- 

log (O/H) 

LOl 

83.9 

5.39 

82.44 ±2.1 

09:54:39.52 

69:05:36.1 

0.54 ±0.03 

0.06 ±0.01 

0.54 

±0.03 

0.74 ±0.03 

1.38 ±0.07 

_ 

_ 

L02 

143.5 

4.73 

-100.43 ±3.3 

09:55:21.58 

69:01:46.9 

0.27 ±0.03 

0.07 ±0.02 

0.57 

±0.03 

0.66 ±0.04 

1.34 ±0.10 

— 

— 

LOS 

59.1 

5.45 

-135.81 ±2.1 

09:55:32.08 

69:01:03.0 

1.32 ±0.05 

0.04 ± 0.02 

0.62 

±0.04 

0.66 ±0.04 

1.33 ±0.11 

-4.23 

-3.54 

L04 

95.3 

6.76 

84.84 ± 1.8 

09:55:34.10 

69:07:29.3 

0.59 ±0.03 

0.11 ±0.01 

0.59 

±0.02 

0.75 ±0.03 

1.44 ±0.06 

— 

— 

L05 

68.8 

7.77 

-199.96 ±0.9 

09:55:47.84 

68:59:28.6 

3.21 ±0.03 

0.18 ±0.01 

0.91 

±0.03 

1.01 ±0.03 

1.16 ±0.04 

-3.83 

-3.08 

L06 

164.9 

12.91 

-217.05 ± 1.5 

09:56:29.29 

68:56:16.5 

2.62 ±0.07 

0.19 ±0.03 

0.53 

±0.03 

0.57 ±0.03 

1.40 ±0.11 

-4.22 

-3.23 

MFOl 

39.7 

5.11 

46.17 ± 1.5 

09:54:44.27 

69:04:23.9 

3.36 ±0.07 

0.11 ±0.03 

0.67 

±0.05 

0.81 ±0.05 

1.42 ± 0.12 

-3.98 

-3.06 

MF04 

29.7 

5.72 

30.58 ± 3.3 

09:54:50.95 

69:02:57.7 

1.49 ± 0.09 

0.12 ±0.05 

0.67 

±0.07 

0.94 ± 0.09 

1.25 ± 0.15 

-4.17 

-3.48 

MF05 

34.0 

9.15 

171.48 ± 2.7 

09:54:54.38 

69:09:19.9 

1.34 ±0.06 

0.22 ±0.03 

0.58 

±0.05 

0.75 ±0.05 

1,44 ±0.14 

-4.26 

-3.55 

MF07 

49.0 

7.06 

159.79 ± 2.7 

09:55:00.02 

69:08:05.5 

1.97 ±0.09 

0.45 ± 0.07 

0.79 

±0.09 

1.14 ±0.12 

1.38 ± 0.17 

-4.04 

-3.34 

MF08 

29.5 

3.72 

104.93 ± 3.3 

09:55:04.24 

69:05:54.6 

2.83 ±0.16 

0.12 ±0.04 

1.06 

±0.19 

1.28 ±0.20 

1.39 ± 0.27 

-3.76 

-3.08 

MFIO 

59.2 

3.63 

-15.29 ± 2.4 

09:55:07.17 

69:03:14.6 

2.34 ±0.06 

0.32 ± 0.03 

0.85 

±0.06 

1.05 ±0.06 

1.34 ±0.10 

-3.99 

-3.27 

MFll 

52.5 

2.44 

39.57 ±3.6 

09:55:09.38 

69:04:14.7 

1.96 ±0.12 

0.38 ±0.07 

0.98 

±0.15 

1.25 ±0.17 

1.43 ± 0.24 

-3.93 

-3.31 

MF12 

40.9 

8.35 

175.98 ± 1.5 

09:55:10.32 

69:08:46.8 

4.84 ± 0.09 

0.12 ±0.03 

0.85 

±0.07 

1.01 ±0.08 

1.42 ± 0.13 

-3.62 

-2.79 

MF16 

57.8 

9.96 

134.91 ± 2.7 

09:55:19.13 

69:09:31.5 

1.97 ±0.09 

0.14 ±0.04 

0.67 

±0.07 

0.94 ± 0.09 

1.44 ±0.17 

-4.09 

-3.32 

MF17 

87.5 

6.30 

166.98 ± 1.5 

09:55:19.72 

69:07:32.6 

1.67 ±0.04 

0.28 ±0.02 

0.96 

±0.04 

1.25 ±0.05 

1.36 ±0.07 

-3.96 

-3.41 

MF19 

55.0 

8.18 

124.71 ± 0.9 

09:55:21.55 

69:08:31.7 

2.16 ±0.03 

0.21 ±0.01 

0.87 

±0.03 

0.98 ±0.03 

1.16 ± 0.04 

-3.98 

-3.29 

MF21 

47.6 

13.46 

-174.18 ± 1.8 

09:55:31.90 

68:56:47.6 

1.36 ±0.04 

0.30 ±0.02 

0.78 

±0.04 

1.00 ±0.05 

1.31 ± 0.08 

-4.12 

-3.55 

MF22 

44.9 

6.37 

-173.58 ±5.7 

09:55:32.27 

69:00:33.3 

0.62 ± 0.09 

0.33 ±0.06 

0.68 

±0.09 

0.95 ±0.11 

1.37 ±0.19 

— 

— 

MF25 

46.6 

6.30 

66.85 ± 1.5 

09:55:42.07 

69:07:00.1 

2.61 ±0.05 

0.22 ±0.02 

0.88 

±0.05 

1.04 ±0.06 

1.25 ± 0.08 

-3.96 

-3.22 

MF26 

56.5 

4.10 

-5.70 ±6.9 

09:55:52.15 

69:05:21.6 

1.09 ±0.13 

0.20 ±0.06 

0.88 

±0.14 

1.16 ±0.17 

1.38 ± 0.24 

-4.04 

-3.56 

MF27 

26.6 

8.04 

-239.23 ± 1.8 

09:55:52.25 

68:59:16.2 

3.04 ±0.07 

0.11 ±0.02 

0.75 

±0.06 

0.88 ±0.06 

1.32 ± 0.11 

-3.96 

-3.12 

MF29 

77.7 

8.33 

-217.95 ±3.9 

09:56:04.49 

68:58:60.0 

0.28 ±0.04 

0.11 ±0.02 

0.49 

±0.03 

0.68 ±0.04 

1.46 ± 0.11 

— 

— 

MF32 

44.0 

5.63 

-191.87 ±3.3 

09:56:15.83 

69:00:51.4 

1.83 ±0.09 

0.46 ± 0.06 

0.85 

±0.09 

1.08 ±0.10 

1.45 ± 0.17 

-4.02 

-3.38 

MF33 

55.3 

4.31 

-100.13 ±2.1 

09:56:16.18 

69:02:39.6 

4.19 ±0.12 

0.22 ±0.05 

0.99 

±0.12 

1.24 ±0.14 

1.30 ±0.18 

-3.59 

-2.84 

MF35 

54.6 

6.58 

-89.64 ±2.7 

09:56:21.77 

69:05:00.9 

0.41 ±0.03 

0.32 ± 0.02 

0.58 

±0.03 

1.09 ±0.05 

1.34 ±0.07 

— 

— 

M82-L1 

— 

69.43 

105.53 ± 1.5 

09:55:58.75 

69:40:09.3 

0.54 ±0.02 

0.05 ±0.01 

0.45 

±0.02 

0.34 ±0.01 

1.29 ± 0.06 

— 

— 

M82-L2 

— 

71.53 

258.72 ± 2.7 

09:55:47.29 

69:41:35.2 

0.55 ±0.02 

0.05 ±0.01 

0.46 

±0.08 

0.30 ±0.05 

1.22 ± 0.18 

— 

— 

M81HII-L1 

— 

23.87 

-169.98 ± 1.5 

09:56:14.29 

68:50:20.7 

1.65 ±0.05 

0.01 ±0.01 

0.19 

±0.01 

0.14 ±0.01 

1.49 ± 0.22 

— 

— 


Abundances derived from the comparison between observation and models on the [N llj/Ho; —[O lilj/Ho; diagram. 
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TABLE 3 

Abundance Gradients of SNRs, H ll regions, and PNe in M81 


Type 

dLog(N/H/)dR 

dLog(0/H)/dR 

References 


SNRs 

-0.023 ± 0.009 

-0.010 ±0.015 

This study 


H II regions 

-0.067 

-0.088 ±0.013 

Stanehellini et al. 

(2014) 


- 

-0.013 ~ -0.020 

Patterson et al. ( 

2oT^ 

PNe 

-0.049 

-0.044 ± 0.007 

Stanehellini et al. 

7^oT4) 


TABLE 4 

Abundance Gradients of SNRs in Nearby Galaxies 


Galaxy 

R 25 [kpc] 

dLog(N/H)/dR 

dLog(0/H)/dR 

References 

M81 

18.34“" 

-0.023 ± 0.009 

-0.010 ±0.015 

This study 

M31 

23.56'" 

-0.040 

±0.004 

Blair & Kirschner (1985) 

M33 

5.77"“ 

-0.089 

-0.140 ±0.040 

-0.035 

-0.090 ± 0.050 

Blair & Kirschner 119851 
Smith et al. (1993) 

MWG 

11.5“' 

-0.09 

- 

Blair & Kirschner (1985) 


^ Ide Vaucouleurs et al.| l |199lD . 

^ Ide Vaucouleurs et al.l il99lD . The distance to M31 adopted is 793 ± 45 kpc JBate et al.l 

f2oa. _ 

IBlair &: Kirschned | |1985D . 
iBlair &: KirschneH dlQSSD . 


TABLE 5 

Oxygen abundance of H ll regions 


12+log(0/H) 


ID 

RGCD[kpc]“" 

PT05' 

KKO 42 

KD023 

B07'* 

M82-Ll'> 

69.43 

8.39 ± 0.04 

9.03 ±0.04 

9.04 ±0.01 

8.58 ±0.01 

M82-L2'> 

71.53 

8.48 ±0.14 

9.06 ±0.17 

9.07 ±0.03 

8.62 ±0.03 

M81HII-L1 

23.87 

8.20 ±0.11 

8.75 ±0.17 

8.73 ±0.02 

8.27 ±0.02 


^ Deprojected distance from the center of M81. 

^ Abundances were derived in the case that M82 objects are fully photo-ionized. 
^ Based on the calibration of IPilvugin fc ThuanI ll2005n . 

^ Based on the calibration of IKobulnickv Sz KewlevI ll2004f) . 

^ Based on the calibration of IKewlev &: Dopital l l2002f) . 

Based on the calibration of IBresolinI <120071) . 



















































